Monopoles, magnetricity and the stray field from spin ice 
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An analysis is presented of the behaviour of muons in the low- temperature state in spin ice. It 
is shown in detail how the behavior observed in some previous muon experiments on spin ice in a 
weak transverse field may result from the macroscopic stray field of magnetized spin ice. A model 
is presented which allows these macroscopic field effects to be simulated and the results agree with 
experiment. The persistent spin dynamics at low temperature originate from the sample and could 
be a muon-induced implantation effect that is operative in out-of-equilibrium systems with long 
relaxation times. 



Spin ice [IH3] is the name given to compounds 
such as Dy2Ti2 07 in which the magnetic Dy ions 
(/i « 10/xb) sit on a pyrochlore lattice (composed 
of corner-sharing tetrahedra) and are constrained by 
easy-axis anisotropy only to point in or out of each 
tetrahedron. The long-ranged dipolar interactions are 
almost perfectly screened at low temperatures [4j [5] so 
that the low-energy properties are essentially identi- 
cal to an effective frustrated nearest-neighbour model 
equivalent to Pauling's model for proton disorder in 
water ice [6]. The effective exchange energy in this 
dipolar Ising magnet is minimized if the moments on 
each tetrahedron satisfy the ice rules, namely that two 
spins point in and two spins point out [lj. This re- 
sults in a state whose large degeneracy is quantita- 
tively consistent with the measured thermodynamic 
properties |7J. Very recently it was realised that the 
natural excitations |8J of spin ice can be described as 
magnetic monopoles |9J : flipping a moment breaks the 
ice rules in two neighbouring tetrahedra (one tetra- 
hedron has moments arranged three-out, one-in; its 
neighbor has one-out, three-in) and these states are 
positive and negative monopoles. Further flips allow 
these monopoles individually to hop through the lat- 
tice. This deconfined monopole picture is an economi- 
cal description of the low-temperature behavior of spin 
ice, as one has moved from considering a concentrated 
collection of localised dipolar spins to a more dilute ar- 
ray of itinerant particles interacting Coulombically, a 
so-called magnetic Coulomb liquid |9J. 

The number of monopoles in spin ice can be esti- 
mated within the framework of Debye-Hiickel (DH) 
theory |10| . see Fig. [TJa), and at very low temper- 
atures the monopoles are very sparse indeed. It has 
been argued [HJ [12] that spin ice should conduct mag- 
netic charge ("magnetricity") in an analogous manner 
to an electrolyte and that Onsager's theory [13J of the 
Wien effect should apply to spin ice. In this context, 
it is useful to consider the key lengthscales: the lattice 
spacing a « 10 A (so that the distance from the cen- 
tre of one tetrahedron to the centre of its neighbour 
is r& = ~ 4.4 A, and the monopole charge is 
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monopoles per tetrahedron) within the framework of 
Debye-Hiickel theory [10]. (b) The key lengthscales in spin 
ice. 



Q = 2/i/r d = 4.50/i B A _1 ), the Bjerrum length £ T = 
fioQ 2 /(87t/cbT) (the distance below which monopoles 
are bound by the Coulomb interaction), the field 
length £b = k^T/QB (the lengthscale above which 
drift of monopoles in a field is discernable) , the mean 
minimum distance between monopoles £ m (related to 
nf, the number density of free monopoles) and the 
Debye length £ B = (k B T/2^ Q 2 n^ 1/2 (the distance 
above which monopole charge density is screened) |12j . 
The temperature dependences of these lengthscales 
are plotted in Fig. [TJb). 

One of the most ingenious experimental demonstra- 
tions of monopoles came from the muon-spin rota- 
tion (/iSR) study of Bramwell et aJ. [llj in which 
Dy2Ti 2 07 was subjected to a weak transverse field 
(TF). At very low temperature (in a regime such that 
£r> ^ £b ^ £t) it was argued that the relaxation rate 
of the muon precession signal A (B) in a magnetic field 
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Figure 2: (a) The extracted monopole charge Q from the 
data in Ref. I 111 showing the plateau close to Q = 5/xbA -1 , 
close to the theoretical value 4.54 Mb A -1 0. The low-T 
data roughly follow a T 2 / 3 dependence, (b) The same data 
transformed to values of m = d ^jf ^ in order to test 
the expected m oc Q 3 /T 2 behavior, showing only limited 
agreement. (c,d) Data replotted from the study of Ref. [14] 
(c) 300 K, filled points for a TF of 7.5 mT, open symbols 
for a longitudinal field of 5mT (which was almost iden- 
tical to zero-field data, but a small longitudinal field was 
used to quench any relaxation due to nuclear spins), (d) 
Measurements in a longitudinal field of 5 mT. Note the re- 
duction in initial asymmetry in the 5K data due to the 
largely static spin ice state. 



Figure 3: (Color online.) Distribution of in the vicin- 
ity of (a-c) two spins at positions (±a,0,0), showing the 
field in the y = (left) and x = (right) planes, and (d-e) 
four spins on the corners of a tetrahedron (inscribed in a 
cube of side length centred at the origin), satisfying the 
ice rules, (d) is shown in the plane z — 0A4d with the 
spins at z = d/2 (z = —d/2) colored black (green), (e) B^ 
in the plane x = 0. The areas with the darkest colors have 
the smallest B^. 



B could be interpreted in terms of a fluctuation rate 
due to monopoles. By analogy with Onsager's theory 
[13J, it was argued that one would expect X(B)/X(0) « 
1 + | + • • • where b = £ T /£ B = VoQ 3 B /8tt/c 2 T 2 . This 
allows one to measure the field dependence of the re- 
laxation rate and then infer the monopole charge Q. 
Bramwell et al. measured an exponentially damped 
precession signal, interpreting the form of the damp- 
ing in terms of magnetic fluctuations resulting from 
monopoles. The analysis involved no free parameters 
and gave a value in very close agreement with the the- 
oretical value Q = 4.50/iBA _1 |9J [see Fig.^a)]. One 
prediction of this method is that the experimentally 
measured quantity m = j(o)^d§^ ~ ^ T~ 2 . As 

shown in Fig. [2^b) , which replots the data of Ref. [TTJ 
m ex T~ 2 only holds for a small subset of the data 
over a very limited range of T. Also, Q oc m 1 / 3 T 2 / 3 
and it is seen in Fig. [2|a) that most of the temper- 
ature dependence in the low-temperature behavior of 
Q is accounted for by the T 2 / 3 factor (present by def- 
inition), so that the closeness of the agreement with 
theory may not be so compelling as it first appears. 
Moreover because the measured m(T) data [Fig.J^b)] 
are found to fall on cooling down to ~ 0.3 K, below 
which they are approximately constant, a limited in- 
tersection with the hypothesized m ex T~ 2 curve is 
not unexpected even if the model is inapplicable. 

Furthermore, despite the elegance of the theoreti- 
cal approach, a surprising feature of these results was 



that any signal in a weak TF could be observed at all. 
An earlier /iSR study using longitudinal-field decou- 
pling [14J had shown that the field at the muon site 
was around 0.5 T, so that a 2 mT field (the maximum 
field used in Ref. |TT]) would not be expected to lead to 
a precession signal. Some data from Ref. 14 are repro- 
duced in Fig. |2jc,d). The relaxation of the precession 
signal in a TF at 300 K is the same as that observed 
in zero field [Fig. [2jc)]. That relaxation rate increases 
dramatically as the sample is cooled (following an ac- 
tivated behavior governed by transitions to/from the 
first excited state doublet of the crystal field [14J). 
At low T [the 5K data are shown in Fig. [2jd)] the 
relaxation is so fast that one can only measure the 
slower relaxation of the remaining "|-tail" (due to the 
component of the muon polarization parallel to the lo- 
cal field). These results are consistent with a largely 
static spin ice state. 

To be sensitive to a weak TF, the field B^ at the 
muon site needs to be very low. The hope that there 
might be some special site of field cancellation is not 
borne out by calculations, and the reason for this 
can be understood from the following considerations. 
First, it is found that the distribution of the magni- 
tude of B^ over the unit cell volume assuming var- 
ious ordered arrangements of spins falls to zero for 
zero B^ [15] . so that if present such zeros are ex- 
tremely isolated and are unlikely to be obtained by 
chance. Even taking just two spins [say at positions 
(±a, 0, 0)] then the field zeros occupy either isolated 
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points [e.g. if the moments \\z, then two zeros are 
found at (0, 0, ±a/y/2), Fig. 3[a)] or in circles [e.g. 
if the moments are along (cos a, 0, ± sin a), there is a 
circle of zeros in the yz plane, centred at (0, 0, § tan a) 

with radius ^2 + | tan 2 a, Fig. 3[b,c)]. For spins on 
the corners of a tetrahedron satisfying the spin ice 
rules only isolated zeros are found [see Fig. J3^d,e)] . 
Crucially their positions do not share the symmetry 
of the tetrahedron but depend on the particular spin 
ice configuration. Since muon sites depend on the 
electrostatic potential and reflect the crystal symme- 
try, then even if a particular muon sits at a site of 
field cancellation, many others will sit at crystallo- 
graphically equivalent sites in which the field is far 
from zero. In a real crystal of Dy2Ti207, containing 
many tetrahedra, the field from neighbouring tetrahe- 
dra (the spins on which can exist in many different 
configurations) will produce additional contributions 
which will displace or remove the zeros. In fact, the 
typical ~ (/xo/47r)/i/r^ ~ 0.1 T, the same order as 
found in the earlier /iSR experiment [14J. This con- 
clusion is in agreement with other recent calculations 
[T6l H7] . Thus it is surprising that applying a TF of 
~ 10 -3 T in the experiment of Ref. [TT] can have had 
any effect at all. Moreover, the recent experiment of 
Ref. [16] provides evidence that no such precession sig- 
nal is in fact observable when Dy2Ti207 is mounted on 
GaAs and the experiment is repeated (so that muons 
missing the sample form muonium and do not con- 
tribute to the observed TF signal). However, the na- 
ture of the signal observed in Ref. [TT| (with the sample 
mounted on silver), which seemingly produces a rea- 
sonable estimate of Q, has remained unaddressed. 

I argue that the most likely resolution of this conun- 
drum is that the TF muon signal which was observed 
originates from outside the sample (as suspected by 
Ref. [TO]) , and furthermore the particular pattern pro- 
duced by a static, macroscopic exterior dipolar field 
produces a muon signal that could be misinterpreted 
as a dynamic signal. Exterior fields are well known 
to result from magnetized ferromagnets but are much 
more unusual in systems with no long range order. A 
macroscopic exterior field directly due to monopoles is 
unlikely because <Csample size (unless T < 0.2 K, 
but even here rif is likely to be larger than predicted by 
DH theory and strongly out of equilibrium [18J, keep- 
ing £d small) but it is more likely due to the spin ice 
magnetization. When this exterior field distribution 
is "sampled" by muons implanted in the silver sam- 
ple holder over the area close to the sample shown 
in Fig. [4^a), simulations show that the resulting re- 
laxation function [Fig. Qb)] mimics the exponentially 
relaxing signal reported in Ref. [TTJ 

The simulation of this phenomenon proceeds as 
follows: In an applied magnetic field B a = /ioH a , 
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Figure 4: (a) The exterior field in the plane of the silver 
plate due to a sample magnetized horizontally and placed 
on top of the plate, evaluated from a closed analytic ex- 
pression. The applied field adds to this exterior field, (b) 
The corresponding muon precession signal evaluated for 
the geometry shown in the inset, together with the 100 mK 
data from Ref. 1111 For this geometry, both x an d the ra- 
tio of sample (grey) to sample holder (pink) were fitted to 
the 2mT data (red open squares). The fitted geometry is 
shown to scale and x — 0.45. The 1 mT line was then eval- 
uated with no further parameter adjustment (data black 
filled circles). 




Figure 5: The field distribution p(B) for the muon aver- 
aged over the silver plate for a variety of sample geometries 
(shown in the insets) and values of %. Also shown are the 
corresponding simulated muon precession signals P z (t) for 
a field of 2 mT. 
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the field inside and outside the sample is given by 
H(r) = (1 — XjV(r))H a , where x is the magnetic sus- 
ceptibility and Af(r) is a symmetrical 3x3 demagneti- 
zation tensor that depends on position r even outside 
the sample, which is of particular interest to us. This 
expression is valid under the assumption that the mag- 
netization inside the sample is uniform. For a sample 
with a cuboidal shape an analytical closed form for 
A/"(r) can be used [19J. With the initial polarization 
along n, the resulting time-dependence of the muon 
polarization can be evaluated using 

(1) 

where the integral is over the area A of the sil- 
ver sample holder surrounding the sample in which 
muons are implanted. Outside the sample B = /i H 
and we have B\\(r) = B n and B±(r) = |B(r) — 
£?ll(r)n|. Thus if B a is along x and one is measur- 
ing P z (t), then B,|(r) = X^ a AT 13 (r) and B ± (r) = 
B a y/ (1 - X-^ii(r)) 2 + X 2 -^ / 2i(r) 2 . The only parame- 
ters in this model are (i) the value of x an d (ii) the 
geometry of the experiment, specifically the size of 
the region outside the sample sampled by the muons. 
Unless x ^> 1, the oscillatory component of P z (t) 
dominates and this can be written J p(B) cos j^Bt dB 
where the function p(B) is plotted in Fig.[5]for various 
geometries and values of x, together with the corre- 
sponding P z {t). The distribution p(B) contains van 
Hove singularities of the macroscopic field distribu- 
tion and broadens with increasing \. It depends on 
the precise shapes and arrangements of crystallites on 
the sample holder, with large density far from B a due 
to regions close to the crystallites. The fit to data 
in Fig. |4jb) is quite reasonable (and corresponds to 
X = 0.45 and the geometry shown in the inset), though 
fits with other geometries are possible. This value of 
X at 0.1K is larger than observed in zero-field cooled 
magnetization data [20J, but falls short of the equi- 
librium susceptibility of spin ice x = 

V37r£ T /r d « 80 

[21J, showing that the magnetization is still slowly de- 
veloping [22] (and this is indeed a manifestation of 
the Wien effect [H]). A single muon measurement 
at a particular B a and T is performed of a timescale 
of ~ 1 hour, so that the effect of slow dynamics (the 
monopole hop rate ~ 10 3 s _1 ) can become important. 
Were larger values of x to be obtained in experiment, 
the simulations in Fig. [5] show that the TF precession 
signals should become strongly distorted by the very 
strong stray fields, an effect that can be looked for in 
future experiments. 

The very fast initial relaxation in the data [llj, not 
fitted in our formalism [see first 1 /is in Fig. |2Jb)] is 
due to spin dynamics from muons stopped in the sam- 
ple |14| . superimposed on the weak TF signal. The 



motion of a monopole is accompanied by the rever- 
sal of a Dy spin, resulting in a change AB^ of up 
to (2/i /47r)2/i/r 3 , where r is the distance from the 
Dy moment to the muon. Within a radius of ~ 50 A 
around the muon there are ~ 10 4 Dy ions, each of 
which can produce a AB^ from tenths to tens of 
mT. Thus even though during the /iSR measurement 
timescale (~ 20 /is) most monopoles appear static 
(and the muon experiences a largely frozen field dis- 
tribution as discussed above), because the monopole 
hop rate is only ~ 10 3 s _1 [20, 23J, the muon is cou- 
pled sufficiently strongly to such a large number of 
Dy ions that even an occasional monopole hop will 
be enough to contribute to measurable relaxation in 
the | -tail. The inferred fluctuation rate of ~ 10 6 s _1 
fT4] from this relaxation is thus entirely consistent 
with this picture. Moreover the muon relaxation rate 
would be expected to track the underlying monople 
dynamics (thus the plateau observed in the relaxation 
time r from a.c. susceptibility [20] is found also in 
/iSR [14J, albeit at a scaled rate). Below 1 K however, 
rif drops dramatically [see Fig. [IJa)] and r diverges, 
while the muon relaxation rate stays approximately 
constant [HI [16 , 24J, an example of the phenomenon of 
persistent spin dynamics observed in many frustrated 
systems [25]. To understand this effect in Dy2Ti2C>7 it 
is worth remembering that the muon itself brings sub- 
stantial kinetic energy at implantation. In the final 
stage of its thermalization, it loses energy from sev- 
eral tens of keV via charge exchange cycles whereby 
an electron successively adds to and is stripped from 
the muon [26J. It is conceivable that this process 
can nucleate monopoles in the sample (the system 
is unable to rapidly transport heat away and is sus- 
ceptible to thermal runaway fZT\ 125]), and when this 
muon-induced concentration exceeds the equilibrium 
concentration, the monopole-induced muon relaxation 
rate will settle at a fixed value. Such an explanation 
could have wider applicability in other frustrated mag- 
nets. 
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